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Possible singlet to triplet pairing transition in Na^.Co02 ■ yH20 
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We present precise measurements of the upper critical field {Hc2) in the recently discovered 
cobalt oxide superconductor. We have found that the critical field has an unusual temperature 
dependence; namely, there is an abrupt change of the slope of Hc2{T) in a weak field regime. In 
order to explain this result we have derived and solved Gor'kov equations on a triangular lattice. 
Our experimental results may be interpreted in terms of the field-induced transition from singlet to 
triplet superconductivity. 



I. INTRODUCTION 

The recent discovery of superconductivity in 
Naa;Co02 • ?/H20)i may provide a unique insight 
into the mechanisms, which determine superconducting 
properties of transition metal oxides. Although the 
superconducting transition temperature, Tc, is much 
lower than Tc's in cuprate superconductors, both system 
share many common features. Co oxide becomes super- 
conducting after hydration that significantly enhances 
the distance between C0O2 layers. This suggests crucial 
importance of the dimensionality. In particular, a quasi 
two-dimensional character of Co oxide shows up in 
the resistivity measurements. Above the transition 
temperature, the in-plane resistivity is three orders 
of magnitude less than out-of-plane one^ Similarly 
to cuprates the Co-based superconductor represents a 
strongly correlated system. The strong correlations may 
be responsible for a nonmonotonic doping dependence 
of Tc- Namely, the critical temperature is maximal for a 
particular carrier concentration and decreases both for 
overdoped and underdoped materials!^ 

However, in contradistinction to cuprates, C0O2 lay- 
ers have a form of a triangular lattice. This feature may 
be responsible for magnetic frustration and unconven- 
tional symmetry of the superconducting order param- 
eter. Investigations of the pairing symmetry with the 
help of nuclear magnetic resonance (NMR)- and nuclear 
quadrupole resonance (NQR)-?^ lead to contradictory re- 
sults. In particular, the presence of nodes in the su- 
perconducting gap remains an open problem. It is also 
unclear whether superconductivity originates from sin- 
glet or triplet pairing. Theoretical investigations do not 
lead to firm conclusions. It has been shown that the 
resonating valence bond state (RVB) may be realized in 
the t-J model on a triangular lattice,^ provided t > Q. 
This may suggest RVB as a straightforward explanation 
of superconductivity in the Co oxides.'' However, it is 
interesting that in addition to singlet superconductivity, 
there is a region of triplet pairing in the phase diagram 
proposed in Refii. Moreover, LDA calculations suggest 
that the ground state of the parent system NaCo204 may 



be ferromagnetic »S Recent density functional calculations 
carried out for Naa;Co204 predict an itinerant ferromag- 
netic state that, however, competes with a weaker anti- 
ferromagnetic instability.^ Triplet superconductivity has 
also been postulated on the basis of symmetry consider- 
ations combined with analysis of experimental results^ 
Therefore, the symmetry of the superconducting or- 
der parameter remains an open problem and both sin- 
glet and triplet pairings should seriously be taken into 
account. In particular, it is possible that singlet and 
triplet types of superconductivity compete with each 
other. In such a case an external magnetic field may 
favor triplet pairing, due to the absence of the Zeeman 
pair breaking mechanism in this state. This should be 
visible in the temperature dependence of the upper criti- 
cal field, Hc2- In order to verify this possibility we carry 
out precise measurements of Hc2- The obtained results 
clearly indicate unconventional temperature dependence 
of Hc2, that cannot be described within the Werthamer- 
Helfand-Hohenberg (WHH) theory. The experimental 
data are compared with theoretical results obtained from 
the solution of the Gor'kov equations on a triangular lat- 
tice. These results may be interpreted in terms of a field- 
induced transition from singlet to triplet superconduc- 
tivity and suggest that phase sensitive measurements to 
distinguish this from other possible interpretations would 
be of great interest. 



II. EXPERIMENTAL RESULTS 

The measurements have been carried out on 
Nao.3Co02 • I.3H2O. Nao.7Co02 (0.5g) was stirred in 
20 ml of a 40x Br2 solution in acetonitrile at room tem- 
perature for 4 days ('Ix' indicates that the amount of Br2 
used is exactly the amount that would theoretically be 
needed to remove all the Na from Nao.7Co02). The prod- 
uct was washed copiously with acetonitrile, followed by 
water and air-dried. After air-drying, the product was 
kept in a sealed container with 100% relative humidity 
for 2 days to obtain the hydrated superconductor. 

All the magnetic measurements were performed using 
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FIG. 1: Temperature dependence of mass magnetization 
M(T) for various magnitudes of magnetic fields. For clar- 
ity of the figure the M{T) curves are offset, except the one 
for OT. The inset presents temperature dependence of the S"^"* 
harmonic susceptibility. 
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FIG. 2: Experimental data for Hc2(T). The horizontal line 
shows the the Clogston-Chandrasekhar limit, whereas the line 
connecting experimental points is only a guide for the eye. 
The insert shows Hc2 (T) in the low field regime fitted by two 
lines. 



DC magnetometer/ AC susceptometer MagLab 2000 Sys- 
tem (Oxford Instruments Ltd.). There is only one super- 
conducting phase transition in the sample and there is no 
coexistence of phases with different critical temperatures 
or critical fields. It has been well confirmed by a single 
peak in a temperature dependence of magnitude of the 
3'''^ harmonic AC susceptibility (see the inset in Fig. 
For third harmonic measurements the AC magnetic field 
of frequency /=lkllz and amplitude ffac — lO^^T was 
applied. For the DC measurements we have applied mag- 
netic fields up to 9T. The temperature was stepped in the 
range about 3-^6K in the case of low and moderate mag- 
netic fields and in the range about 2-^5K for the highest 
fields. The size of the step was 20mK and the temper- 
ature was stabilized during each measurement with the 
accuracy 2mK. A set of typical M{T) curves recorded for 
applied fields of 0, 2, 4, 6, and 8T is shown in Fig. ^ 
Except the case of lowest magnetic fields the magneti- 
zation is positive in the whole temperature range. It is 
due to the domination of ferromagnetic and/or paramag- 
netic contributions in the total magnetization at higher 
fields. The superconducting transition manifests itself as 
a downturn in M(T) at low magnetic fields, whereas at 
the higher fields only the change in the slope in M{T) 
is observed. This enabled a simple determination of the 
critical temperature; namely, Tc was determined from the 
intersection of the two straight -lines that fit relevant lin- 
ear regimes (see Fig. The zero-field critical temper- 
ature determined in this way is Tc(0) = 4.345 ± 0.015 K. 
The results of the measurements are presented in Fig. |5] 
in equivalent form as Hc2{T). 

Close to Tc{H = 0) one can expect that the Ginzburg- 



Landau theory gives accurate results and, therefore, tem- 
perature dependence of Hc2 should be linear. In cuprate 
superconductors Hc2 shows unconventional temperature 
dependences^ Close to TdH = 0) 7?c2(r) is almost 
linear and, then, the curvature smoothly increases with 
the decreasing temperature. However, as can be inferred 
from Fig. [3 this is not the case for Naa;Co02 -2/1120. For 
IT < _ff < 3T the experimental data can be fitted very 
well by a linear function. However, such a fit deviates 
from experimental points for weaker magnetic field. Sim- 
ilar temperature dependence of Hc2 has been obtained, 
e.g., from the specific heat measurements For H < 3T 
the experimental data presented in Fig. |21 as well as 
those reported in Ref>i^, can be fitted by two linear 
functions. In our case they are: Hc2{T) = 7.4 — 1.7T 
in the weak field regime and Hc2(T) = 40 — 9.4T for 
stronger magnetic field. Using the WHH formulajii 
Hc2{0) — 0.7Tc {dHc2/dT) \t^, one estimates correspond- 
ing values of i?c2(0)'s equal to 5.2T and 28T, respectively. 
Such a behavior may originate from competition between 
two superconducting order parameters with close transi- 
tion temperatures but different temperature dependences 
of Hc2- Singlet and triplet order parameters are possible 
candidates due to the absence of Zeeman pair breaking 
in the latter case. i?c2(0) obtained from the WHH for- 
mula in the strong field regime is beyond the Clogston- 
Chandrasekhar (CC) limit. i7c2(0)'s reported in Ref.A^ 
and estimated from Reim^ are even higher. Although, 
the extrapolated 77c2(0) may be overestimated, our ex- 
perimental data clearly show that Hc2 exceeds the CC 
limit aheady for T ~ 0.6 T^. The large slope of Hc2{T) 
suggests that even in the case of renormalization of the 



paramagnetic pair breaking mechanism (e.g., similar to 
that in the strong-coupHng electron-phonon approacbiSi) 
Hc2{0) should be beyond the CC limit. This speaks in 
favor of triplet superconductivity. On the other hand, 
Hc2{0) estimated from the low field data does not ex- 
ceed CC limit. Therefore, superconductivity in a weak 
magnetic field may originate from the singlet pairing. In 
the following we show that this tempting interpretation 
of experimental data remains in agreement with theoret- 
ical results obtained from the numerical solution of the 
Gor'kov equations. Our fit neglects a positive curvature 
of Hc2{T) that occurs for H < 0.9 T. At the end of this 
paper, we discuss possible origins of this feature. 
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FIG. 3: Order parameters for superconductivity on a trian- 
gular lattice. Such a set of order parameters is introduced in 
each pairing channel, i.e., for A, A^^, A^^, A-^-^. 



III. THEORETICAL APPROACH TO THE 
UPPER CRITICAL FIELD 

In order to calculate the upper critical field we consider 
a triangular lattice immersed in a uniform perpendicular 
magnetic field: 

i,cr 

i 

+ E E (A--cL^ct^^+h.c.). (1) 

(ij) (Ti,cr2 = Ti 

tij is the hopping integral between the sites i and j in 
the absence of magnetic field and 6'y is the Peierls phase 
factor, responsible for the diamagnetic response of the 
system: 

% = — / A-dl, (2) 

where $o = hc/e is the flux quantum. The chemical 
potential /x has been introduced in order to control the 
carrier concentration. In the Hamiltonian Q 

Ay = (c,|Cjx - c^iCji) (3) 

and 

AI/ = (c^^cji+c^icj^), (4) 

Ai; - (qtc,t>, (5) 

a!/ = (cac.i) (6) 

denote the pairing amplitudes in singlet and triplet chan- 
nels, respectively. 

In order to determine the upper critical field we follow 
the procedure introduced in Refs»i2*i^. Namely, we diag- 
onalize the kinetic part of the Hamiltonian [the first term 



in Eq. (Q] by introducing a new set of fcrmionic opera- 
tors. In the Landau gauge A = (—i/^y, 0,0) this set of 
fermionic operators is determined by a one-dimensional 
eigenproblem known as the Harper equation. We restrict 
further considerations to the nearest-neighbor hopping, 
i.e., Uj — t for the neighboring sites i,j and otherwise. 

Using the solutions of the Harper equation we write 
down a self-consistent equation for the gap functions. 
Magnetic field breaks the translational symmetry and, 
therefore, the order parameters are site dependent. How- 
ever, in the chosen gauge they depend on the y coor- 
dinate only. As the superconductivity develops on the 
triangular lattice, at each site we introduce three order 
parameters (A^, A^, A^) in each pairing channel, i.e., for 
A, A^-'-, A^1^, A^J- (see Fig. O. In order not to assume 
any particular pairing symmetry we consider these order 
parameters as independent quantities. In the following 
we do not assume any particular orbital symmetry of the 
pair state. However, independently of this symmetry all 
these order parameters vanish at Tc- Therefore, the gap 
equation can be expressed with the help of three vec- 
tors Ai'2.3^ where A" = {A^, A^, . . .). The lower index 
enumerates rows of the lattice sites, whereas the upper 
one indicates the direction of the bond, as depicted in 
Fig. |31 Hc2 is defined as a field, at which all components 
of these vectors vanish. This can be determined from 
the linearized version of the gap equation that is of the 
following form: 




For the sake of brevity we do not present an explicit 
form of A4. This matrix can be expressed with the help 
of the Cooper pair susceptibility and eigenfunctions ob- 
tained from the Harper equation. 

The temperature dependence of Hc2 has been obtained 
from Eq. Q for singlet and triplet superconductivity. In 
the latter case we have investigated separately the paired 
states lii), 1/V2{\U)+ lit)), and |TT). We refer to these 
states by the corresponding spin projection 5'z = — 1, 0, 1, 
respectively. These states are affected by the magnetic 
field in different ways due to the Zeeman coupling. In 
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FIG. 4: Numerical results for hc2{T). Presented results have 
been obtained for occupation number n — 0.95 and for V" = 
0.55t, V* = 0.75t. 
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FIG. 5: Numerical results for hc2{T). We have assumed the 
model parameters n = 0.67 and V" = V* = 0.7t. For these 
model parameters kTc — 0.2t. The experimental points have 
been taken from Ref^ 



the case of equal-spin-pairing this coupling is almost in- 
effective as it leads to a renornialization of the chemical 
potential only /i^/i = /i±l/2 g/isiJ^. 

Figures 01 and O show the numerical results obtained 
for 150 X 150 cluster. In particular, in Fig. |S1 we have 
chosen the occupation number that is close to the ex- 
perimentally determined optimal doping. We present a 
reduced magnetic field h = 27r<f>/<i>o, where $o is the flux 
quantum and $ is the magnetic flux through the lattice 
cell. When the zero-field transition temperatures for sin- 



glet and triplet pairings are of comparable magnitudes, 
the slope of Hc2{T) is much larger for triplet supercon- 
ductivity. As a result, the triplet superconductivity is 
characterized by much larger Hc2{T — 0). It is interest- 
ing that this feature remains valid also for triplet pairing 
with Sz ^ that is affected by the Zeeman pair break- 
ing (see Fig. As expected, the highest value of the 
upper critical field is obtained for the triplet equal-spin- 
pairing. Comparing Figs. 0] and [Hone can see that the 
above statements on the slope of Hc2(T) hold in a wide 
range of the model parameters. However, due to limi- 
tation of the cluster approach it is difficult to perform 
numerical calculations at very low temperatures^ From 
Figs. O and 121 follows that experimental data are qualita- 
tively reproduced when Tc{H — 0) for triplet supercon- 
ductivity is slightly less than the transition temperature 
for the singlet one. Then, sufficiently strong magnetic 
field leads to a transition from singlet to triplet super- 
conductivity that shows up in a change of the slope of 
Hc2{T). 

The external magnetic field affects the relative phase 
of the order parameter in different directions presented in 
Fig. O According to Eq. Q this phase can change from 
site to site and, therefore, it is impossible to determine 
globally the type of the symmetry of the energy gap. 
However, we have found that for singlet pairing Tc{H = 
0) is exactly the same as transition temperature obtained 
for di + id2 symmetry, according to the notation in Ref 
On the other hand, for the triplet pairing Tc{H — 0) 
corresponds to that for the /-wave pairing when ri < 1 
and px + ipy symmetry for larger occupation number. 
The exact position of the boundary between both the 
triplet solutions depends on the pairing potential. 



IV. DISCUSSION AND CONCLUSIONS 

Our linear fit to the experimental data in the inter- 
mediate field regime seems to be very accurate, strongly 
supporting the triplet pairing. Even stronger evidence 
comes from the presence of superconductivity above the 
CC limit. However, at weak field there exists also other 
possibility: Hc2{T) for the field less then approx. 0.9T 
can be fitted by a concave curve. Similarity between 
Naa;Co02 • 1/H20 and high-Tj. superconductors may sug- 
gest a common mechanism that leads to the positive cur- 
vature of Hc2{T). This also speaks in favor of a sin- 
glet pairing in this regime. In such an approach the 
singlet-triplet transition takes place at slightly higher 
field, H w 0.9 T. From a theoretical point of view the 
upward curvature of Hc2(T) can occur for instance in: ex- 
tremely type II superconductors described by the boson- 
fermion model, the systems with a strong disorder suf- 
ficiently close to the metal-insulator transitionjSS the dis- 
ordered superconductors due to mesoscopic fluctuation^Si 
Josephson tunneling between superconducting clusterS)SS 
in a mean-field-type theory of Hc2 with a strong spin-flip 
scattering, and due to a reduction of the diamagnetic 
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pair-breaking in the stripe phase.--- Other theoretical ap- 
proaches to this problem include, e.g., the superconduc- 
tivity with a mixed symmetry (s + d) order parameteB^^ 
and Bose-Einstein condensation of charged bosonsi^ 

Hc2{T) obtained from the resistivity measurements^l 
is lower than presented here and, e.g., in Refsii^^. In 
particular, it is lower than Hc2{T) obtained from mag- 
netization measurements on quasi-single crystals. This 
discrepancy remains unexplained. One of possible expla- 
nations is that it originates from the presence of lattice 
defects, that for short coherence length superconductors 
form Josephson junctions. These junctions affect the re- 
sistivity measurements much stronger than the magneti- 
zation ones. 

To conclude, we have measured the temper- 
ature dependence of the upper critical field in 
Nao.3Co02 • I.3H2O and we have found an interesting 
feature; namely, an abrupt change of slope of Hc2{T) in 
a weak-field regime. This feature is in qualitative agree- 
ment with results reported in Refii^. Moreover, such a 
bend in a weak field regime is visible also in other magne- 
tization measurements in specific heat measurements^ 
as well as in resistivity measurements (see Fig. 4a in 
Ref^.) In order to explain the origin of such a behavior 
we have solved Gor'kov equations on a triangular lattice 
for singlet and triplet types of pairing. Our experimen- 



tal results are consistent with a scenario of competing 
singlet and triplet superconductivity. Within such an 
approach magnetic field induces a transition from singlet 
to triplet superconductivity that shows up in a change 
of slope of Hc2{T): in a weak magnetic field the sin- 
glet pairing takes place and sufficiently strong magnetic 
field drives the system into the triplet state. Recently, 
a field-induced transition between various types of sin- 
glet superconductivity has been proposed to take place 
in cuprates^ (an occurrence of a minor idxy component 
of the order parameter). 
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